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摘 ， 要 : 碳 源 / 汇 是 解释 地 球 大 气 碳 循环 过 程 的 重要 指标 ,探究 三 江 源 的 碳 源 / 汇 特征 对 于 理解 该 地 
区 植被 对 全 球 气候 变化 的 响应 具有 重要 意义。 三 江 源 以 脆弱 的 草地 生态 系统 为 主 , 且 对 全 球 气候 
变化 非常 敏感 。 该 地 区 生态 环境 极其 脆弱 ,大 部 分 地 区 条 件 恶 劣 导致 实测 数据 稀缺 ,很 难 对 该 地 
区 的 碳 源 / 汇 时 空格 局 进行 完整 剖析 。 因 此 通过 以 三 江 源 5 种 典型 草地 群落 (人 金 露 梅 紫花 针 茅 、 风 
毛 菊 、 小 某 草 、 及 青藏 草草 群落 ) 为 研究 对 象 ,基于 BIOME-BGC 模型 ,利用 地 理 数 据 、 气 象 数据 和 植 
被 生理 参数 等 数据 ,得 出 2001 一 2017 年 三 江 源 草地 群落 的 净 初 级 生产 力 (NPP) 净 生 态 系统 生产 
力 (NEP) 模 拟 值 ,并 对 草地 群落 NPP NEP 交 化 特征 与 气温 、 降 水 相关 性 以 及 碳 利用 效率 变化 等 特 
征 进行 了 综合 分 析 。 结 果 表 明 : 三 江 源 区 NPP、NEP 在 空间 格局 上 ,表现 为 由 东南 向 西北 数值 逐渐 
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中 ,人 金 露 梅 群落 NPP 平 均值 最 高 为 342.00 g Cem .a ,青藏 草草 群落 NPP 平 均值 最 低 为 55.93 gC 
m +a "1;5 种 草地 群落 NEP 的 多 年 平均 值 为 49.02gC.m”.a , SEM KES RFR SSI 
被 群落 的 WEP 值 呈 缓 慢 的 上 升 趋势 , 风 毛 菊 和 小 萝 草 群落 呈 缓 慢 下 降 趋 势 。 研 究 发 现 三 江 源 草地 
生态 系统 具有 显著 的 碳 汇 作用 , 且 不 同 群 落 NPP、VEP 对 气温 和 降水 的 响应 程度 有 所 差异 ,5 种 群 
落 NPP 与 气温 均 哇 显著 正 相关 ,但 NPP、NEP 与 降水 量 的 相关 性 较 低 ;5 种 群落 均 具 有 和 较 强 固 碳 潜 
A , 除 金 露 梅 外 其 余 植 被 群落 的 碳 利用 率 均 在 0.625 以 上 。 

关 键 词 : 碳 源 / 汇 ; NPP; NEP; BIOME-BGC 模型， 气候 因子 ; 三 江 源 


文章 编号 : 


草地 作为 自然 生态 系统 不 可 或 缺 的 一 部 分 , 它 
有 防风 固沙 、 养 育 生物 、 涵 养 水 源 等 生态 功能 ,不 
仅 能 够 促进 生态 圈 的 物质 和 能 量 循环 ,还 具有 固 碳 
能 力 ,能 够 调节 全 球 碳 平衡 .维持 全 球 气候 稳定 六。 
随 着 全 球 社会 经 济 快速 发 展 ,石油 .煤矿 等 自然 资 
源 快速 消耗 ,地 球 表 面 CO; 浓 度 也 不 断 攀 升 。 据 相 
关 数 据 表 明 ,2013 4E CO, HR BE (7.20 10 kg +m?) ït 
高 于 工业 革命 前 的 浓度 (5.04x107 kgm) ,被 称 为 
全 球 变 暧 “新 时 代 ” 的 起 点 。 自 工业 化 以 来 ,全 球 的 
平均 气温 上 升 了 约 1% ,温室 效应 日 渐 严 重 , 这 一 系 
列 的 问题 都 受到 国内 外 学 者 的 广泛 关注 。 净 初级 


= 
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生产 力 (Net Primary Productivity , VPP) 又 名 净 第 一 
性 生产 力 , 是 植物 通过 光合 作用 吸收 CO, .固定 并 转 
化 为 有 机 碳 所 形成 的 能 量 减 去 自身 用 于 异 氧 呼吸 
的 能 量 , 能 够 反映 植物 的 生长 状况 .区域 生态 系统 
的 健康 状况 ,其 至 是 表征 全 球 陆 地 碳 平衡 和 碳 循 环 
过 程 的 重要 指标 。 净 生态 系统 生产 力 (Net Eco- 
system Productivity , NEP) 则 是 NPP 减 去 异 养 呼吸 的 
部 分 ,反映 生态 系统 的 物质 循环 以 及 能 量 流动 ,其 
数值 大 小 直接 反映 碳 源 / 汇 的 大 小 ,是 表征 碳 源 / 汇 
的 重要 指标 ”。 碳 源 即 植物 向 外 界 排放 CO,, 碳 汇 即 
植物 吸收 空气 中 的 CO， 。 
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因此 , 自 20 世 纪 60 年 代 “ 国 际 生 物 圈 计划 ” 实 
施 以 来 ,专家 学 者 对 全 球 碳 循环 过 程 中 的 关键 指标 
以 及 影响 生态 系统 碳 循环 的 因子 进行 了 大 量 测算 
和 人 研究。 测算 植被 生物 量 的 方法 主要 有 模型 模拟 
法 和 传统 实测 法 。 传 统 实测 法 在 研究 初期 运用 较 
多 ,是 基于 研究 区 域 的 实测 数据 来 估算 草地 生物 
量 ,方法 简单 .精度 高 ,但 具有 破坏 性 , 且 只 适用 于 
小 面积 研究 区 ,在 条 件 恶劣 人迹 罕 至 的 地 区 该 方 
法 较 难 实施 “。 因 此 ,利用 模型 估算 NPP NEP WURA 
一 种 广泛 应 用 的 重要 手段 。 模 型 模拟 方法 随后 得 到 
快速 扩展 ,最 初 多 应 用 统计 模型 ,然后 发 展 为 结合 遥 
感 数据 的 光 能 利用 模型 .生态 过 程 模 型 ,如 MAIAI7、 
CASA”! YPM"! BIOME-BGC2 BEPS2 等 模型 。 
总 之 , 碳 循 环 过 程 的 研究 正在 不 断 扩 充 ,研究 方法 
和 手段 亦 在 不 断 完善 ,从 传统 实地 测量 到 多 技术 应 
用 ,从 单一 实测 数据 到 多 源 数据 融合 ,研究 对 象 也 
从 局 部 小 范围 区 域 上 升 至 全 球 区 域 '”。 

从 我 国 近 千年 的 历史 资料 来 看 ,青藏 高 原 的 气 
候 变 化 比 我 国 东部 的 变化 要 早 10 ~ 60 a ,被 称 为 
气候 变化 的 启动 区 ,其 气候 变化 的 幅度 也 远 高 于 其 
它 地 区 ,因此 是 全 球 气候 变化 的 敏感 区 之 一 "”。 位 
于 青藏 高 原 腹 地 的 三 江 源 为 寒冷 干旱 生境 ,其 植被 
处 于 早期 发 育 期 ,系统 功能 结构 较为 简单 ,以 草地 
生态 系统 为 主 ,生态 环境 十 分 脆弱 ,一旦 遭 到 破坏 
将 很 难 恢 复 。 而 关于 三 江 源 生态 环境 的 理论 和 实 
践 知 识 还 较为 缺乏 , 且 随 着 人 类 足迹 在 三 江 源 的 增 


多 ,势必 会 给 该 地 区 的 生态 环境 带 来 一 系列 的 问 
题 。 基 于 此 ,本 文 在 三 江 源 区 选取 了 5 种 典型 草地 
群落 ,通过 BIOME-BGC 模 型 ,得 到 了 2001 一 2017 年 
草地 植被 的 模拟 NPP NEP 值 ;并 分 析 了 5 种 草地 和 群 
落 的 碳 / 源 汇 动态 变化 特征 、 植 被 碳 利 用 效率 ,并 评 
佑 气候 变化 对 草地 碳 / 源 汇 的 影响 ,揭示 了 三 江 源 草 
地 生态 群落 的 固 碳 潜力 对 气候 变化 的 响应 机 制 , 为 


客观 评价 三 江 源 生态 系统 的 碳 源 / 汇 功 能 和 开展 实 
践 活动 提供 参考 。 


1 研究 区 域 .数据 方法 


1.1 研究 区 域 

三 江 源 地 处 青海 省 南部 “世界 屋 疹 "青藏 高 原 的 
腹地 ,地 理 位 置 在 31° 30’ ~36° 12' N, 89°45’ ~ 
102 "23 外 之 间 ,是 长 江 .黄河 以 及 澜沧江 的 发 源 
地 。 主 要 地 貌 类 型 为 山地 ,平均 海拔 为 3 500~4 
800 m。 三 江 源 总 土地 面积 为 30.25x10' km ,其 中 
草地 面积 占 全 区 总 面积 的 65.37%。 不同 地 区 干 
湿 差 异 显 著 "“ , 随 着 地 域 由 东南 向 西北 延伸 ,海拔 
逐渐 升 高 .水 热 组 合 也 逐渐 降低 ,该 地 区 的 气温 和 
降水 历年 平均 值 分 别 为 -5.38~4.14 C 、262.2 ~ 
772.8 mm, 具 有 独特 的 高 寒 生 态 系 统 。 根 据 三 江 源 
的 植被 分 布 (图 1) ,本 文选 取 了 5 种 典型 草地 植被 及 
其 所 对 应 的 7 个 站 点 ,依次 为 : 金 露 梅 群落 ( 班 玛 
H) VNE RIE OA A H .河南 县 ) 风 毛 菊 群 落 ( 杂 


0 200 km 
CS， 


o 选择 站 点 m 高 山 植被 ” 咖 阔 叶 林 、 针 叶 林 
en ee m 其 他 植被 
= 


无 植被 地 段 
m 灌 从 


图 1 三 江 源 地 区 植被 及 站 点 分 布 示意 图 


Fig. 1 Location of vegetation and meteorological site in the Three River Source Area 
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多 县 ) 紫花 针 茅 群落 ( 曲 麻 莱 县 .五 道 染 ) 和 青藏 蔓 
草 群落 ( 沱 沱 河 )。 
1.2 数据 来 源 

研究 所 需 植被 数据 来 源 于 1980 年 1:100 万 全 国 
植被 分 类 数据 ,包括 植被 型 .植被 型 组 .植被 群 系 和 
植被 大 类 ,最 终 确 定 研 究 对 象 为 5 种 草地 植被 。 站 
点 数据 来 源 于 全 国 气 象 站 点 基础 数据 ,根据 研究 区 
草地 分 布 状况 ,选取 研究 区 及 其 周边 的 7 个 观测 站 
点 及 4 个 基准 站 点 ( 达 日 GRE ETE .格尔木 ) ,主要 
提取 它们 的 高 程 ` 经 纬度 信息 。 气 象 数据 来 源 于 中 
国 气象 数据 网 站 (http://ede.cma.gov.cn) ,下载 各 个 站 
点 2001 一 2017 年 的 全 国 地 表 数 据 , 包 含 日 最 高 温 
度 .日 最 低温 度 .日 平均 温度 日 总 降水 量 。 利 用 山 
地 气候 模拟 器 (MT-CLIM ) 模 拟 得 到 站 点 的 饱和 水 汽 
压 亏 缺 值 .地 表 太 阳 辆 射 .日 照 长 度 , 作 为 BIOME- 
BGC 模 型 的 输入 参数 ,植被 生理 参数 数据 通过 查阅 
文献 获得 , 表 1 为 本 文中 模型 的 主要 参数 设置 。 
1.3 BIOME-BGC 模型 

Biome-BGC 模型 是 由 FOREST-BGC 模型 不 断 完 
善 而 形成 的 一 个 生物 地 球 化 学 模型 。 该 模型 基于 
气象 数据 (气温 .降水 .C0, 浓 度 ) .基础 地 理 数据 以 
及 植被 生理 参数 数据 ,利用 植被 的 光合 作用 、 自 养 
呼吸 . 异 养 呼吸 .土壤 凋落 物 等 模型 的 原理 ,以 天 为 
单位 ,通过 模拟 碳水 、 氮 这 3 个 关键 因子 的 循环 , 进 
而 模拟 植被 每 天 的 生长 状况 。 因 此 ,根据 模型 输出 
参数 可 以 估算 植被 乃至 生态 系统 的 生产 力 。 模 型 


表 1 植被 生理 生态 参数 表 


Tab.1 Physiological parameters of vegetation 


参数 C3 草地 参考 文献 
FEW HERE H 1 
凋落 期 占 生 长 季 比 侦 1 
年 叶 和 细 根 的 转换 比例 1 
落叶 层 各 部 分 的 比 侦 0.39 
叶 凋 落 物 的 纤维 素 比 例 0.44 
叶 凋 落 物 木质 素 比 俩 0.17 
细 根 各 部 分 的 比例 0.30 WHITE" 
细 根 的 纤维 素 比 例 0.45 WHITE 
细 根 木质 素 比 例 0.25 WHITE" 
枯木 纤维 素 比 例 0.75 REEVES"! 
枯木 木质 素 比例 0.25 REEVES" 
水 汽 压 差 限制 传导 上 限 930 
水 汽 压 差 限制 传导 下 限 4100 


的 运行 数据 和 参数 包括 3 部 分 :植被 生理 生态 参数 
(epc) 初始 化 文件 (ini) .日 值 气象 数据 (metdata ) 。 
模型 运行 的 主要 公式 ”为 : 
(1) 气 孔 导 度 ( g, ) 的 计算 : 
E, = Gna X SE(T an)” SE(VPD) X SF) (1) 
式 中 : g,, 为 最 大 气孔 导 度 ; SF 为 极限 应 力 函 数 ; 
VPD 为 水 汽 压 亏 缺 ; 7 为 最 低温 度 ; y HERR 
水 量 。 
(2) 光 合作 用 速率 (4) 计 算 : 
A=min(A,,A)) (2) 
式 中 : 4, 为 RuBP( 植 物 光 和 作用 中 起 重要 作用 的 一 
种 五 碳 糖 ) 饱 和 时 的 光合 作用 效率 ,由 学 者 
(Faquthar , VonCaemmerer) 研 究 建 立 : 
Vanai Ty) 


C, efi + z) 
式 中 :Vi 为 Rubisco 送 化 酶 (光合 作用 中 C3 反应 
时 起 重要 作用 的 一 种 羧 化 酶 ) 的 最 大 羧 化 能 力 ; C; 
为 叶片 内 的 CO; 浓度; 为 CO: 的 光 补 偿 点 ; 天. 为 
C0; 的 米 氏 常数 ; p 为 植物 叶片 内 0; 分 压 ; R, 为 在 
光照 下 叶片 的 呼吸 速率 ; k, 为 针对 CO: 的 竞争 抑制 
常数 。 
最 大 法 化 能 力 (V,,, ) 的 计算 如 下 : 


P act X flnr 
cmax fnr x sla x C:Nleaf 


式 中 : act X Rubisco ff; flnr, far 分 别 为 叶 氮 BF 
Fr EATE Rubisco 酶 分 子 中 所 分 配 的 比例 ; sla 为 实 
际 的 叶片 面积 ; C:Nleaf 为 叶片 中 的 碳 与 氮 的 总 含 
量 之 比 。 
公式 (2) 中 4 为 辆 射 限制 下 的 光合 作用 速率 ; 
J(C -TD 


4 = -R, (3) 


(4) 


i 95,4 1050, a 
式 中 :为 RuBP 再 生 速 率 。 
(3) 自 养 呼吸 速率 (R) 计 算 : 
1-20 
R=R Qo" , R.=leafn X Mrpern, (6) 


Mrpern = 0.218(kgC/kgN /a) 


式 中 ;Mrpern JIE AA Œ ; leafn JIU Hr AY BUR E ; 
R, 为 环境 温度 在 20 % 时 植被 的 呼吸 速率 。 
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2 结果 与 分 析 


2.1 模型 结果 验证 

表 2 为 模拟 NPP 数 据 与 实测 数据 的 比较 ,由 于 
实测 点 与 模拟 站 点 并 不 完全 重合 ,上 且 有 的 模拟 站 点 
附近 条 件 较 为 恶劣 ,因此 缺少 实测 值 。 根 据 表格 中 
NPP 值 的 比较 可 以 发 现 :2017 年 56046 站 点 附近 的 
实测 值 远 高 于 模拟 值 ,原因 是 此 站 点 位 于 巴 颜 喀 拉 
山脉 ,海拔 高 于 最 近 的 实测 点 。 同 时 2017 年 560018 
站 点 模拟 值 与 实测 点 之 间距 离 较 远 ,导致 实测 值 低 
于 模拟 值 ; 除 此 之 外 ,其 它 几 组 数值 相差 较 小 。 

通过 模拟 NEP 数 据 与 BEPS 模型 模拟 值 比 较 
( 表 3) ,验证 时 期 均 为 7.8 月 份 的 数据 。 其 中 56046 
站 点 与 验证 数据 相差 较 小 ,其 余 3 个 站 点 值 均 略 高 
于 验证 值 。 说明 本 模型 的 模拟 值 与 真实 值 较为 接 
近 ,模拟 的 精度 较 高 ,能 够 对 三 江 源 的 碳 源 / 汇 进行 
人 研究 。 
2.2 NPP、 VEP 趋 势 及 空间 特征 分 析 
2.2.1 NPP NEP 日 变化 趋势 ”三 江 源 5 种 典型 草 
地 植被 NPP NEP 的 时 间 变 化 特征 (以 2017 年 数据 
为 代表 ) 表 现 出 由 东南 向 西北 走向 (图 2a ~ j) ,植被 
生长 活跃 期 逐渐 缩短 ,由 250 qd 减少 到 130 d。 在 年 
际 上 ,NPP NEP 数值 均 表现 为 : “春季 平 稳 -夏季 大 
幅 上 升 - 秋 季 大 幅 下 降 - 冬 季 平 稳 ” 的 倒 “U” 型 趋势 ， 


表征 出 “ 碳 源 - 碳 汇 - 碳 源 ” 特 征 。 日 NPP 峰 值 由 高 
到 低 排 序 依次 为 : 金 露 梅 (2.1$ g C+-m a >N 
(1.85 ¢C+m™?+a')> UE (1.61 gC:m”…a"')> 紫 花 针 
A (1.19 gCom?+a")> Fy CE AL (0.8 gCem™ a") & 
春 两 季 植 被 NEP 为 负 值 ,主要 是 因为 此 段 时 期 植被 
的 生长 呼吸 作用 大 于 植被 的 光合 作用 ;开春 后 气温 
回升 ,水 热 供 给 增多 植被 开始 进入 生长 期 ,NEP 剧 
烈 上 升 ,秋季 叶片 开始 凋落 ,植物 开始 进入 休眠 期 ， 
NEP 则 剧烈 下 降 。 

2.2.2 NPP NEP $ REM 三 江 源 5 种 草地 和 群落 
2001 一 2017 年 WPP 均 值 变化 得 出 (图 3a ~ e) ,其 变 
化 均 维 持 稳定 上 升 趋势 ,年 增长 速率 依次 为 金 露 梅 
(5.63% )> 小 蒿 草 (5.21% )> 风 毛 菊 (4.84% ) > 紫花 针 
F (4.55% ) > 青藏 草草 (1.12% ) 。 较 为 严寒 的 地 区 增 
长 速率 较 低 ,而 气候 条 件 较 好 的 地 区 增长 率 较 大 ， 
其 中 青藏 蔓草 NPP 年 平均 值 上 升 缓慢 ,处 于 47 ~ 
68 gC:m ”a 之 间 , 原 因 是 该 地 处 于 昆仑 山脉 , 气 
候 寒 冷 ,气温 对 于 植被 生长 的 影响 大 于 降水 的 影 
响 , 即 青藏 蔓草 的 生长 对 气温 的 响应 程度 要 高 于 对 
降水 的 响应 程度 ,气温 是 植被 生长 的 限制 因子 2.1。 

5 种 群落 17 a 来 的 NEP 值 分 析 得 出 (图 3f ~j))， 
金霞 梅 、 紫 花 针 茅 和 青藏 莹 草 的 NEP 值 增长 幅度 较 
小 ,小 葛 草 和 风 毛 菊 总 体 呈 缓慢 下 降 趋 势 。5 种 草 
地 群落 的 WEP 值 均 大 于 0, 可 见 三 江 源 区 的 草地 和 群 


表 2 净 初 级 生产 力 (NPP) 的 模拟 值 与 实测 值 的 验证 比较 


Tab.2 Verification of simulated NPP with observed values 


年 从 模拟 数据 实测 数据 
站 点 编号 经 度 /?E 纬度 /*N 模拟 NPP 经 度 /?E 纬度 /"N 实测 NPP 
2017 56046 99.65 33.15 225.119 99.84 33.60 435.118 
2017 56065 101.60 34.73 222.396 100.92 35.32 222.075 
2017 56021 95.80 34.12 156.578 95.81 34.02 123.795 
2017 56018 95.28 32.88 218.537 95.56 33.21 102.758 
2017 56151 100.75 32.93 323.684 100.57 35.78 395.483 
2012 56004 92.43 34.22 50.129 92.49 34.27 65.205 
2005 56004 92.43 34.22 43.131 92.43 34.28 65.227 
表 3 净 生 态 系统 生产 力 (NEP) 的 模拟 值 与 BEPS 模 型 模拟 值 的 验证 比较 

Tab.3 Validation and comparison of NEP simulation value and BEPS model simulation value 
年 份 站 点 编号 ARENE 纬度 /eN 本 文 模拟 NEP 值 BEPS 模 型 模拟 NEP 值 
2004 56004 92.43 34.22 15.646 2.485 
2004 56018 95.28 32.88 35.742 18.611 
2004 56021 95.80 34.12 32.904 22.734 
2004 56046 99.65 33.75 38.591 39.189 
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(a) 金 露 梅 VPP (b) 小 草草 NPP (c) 风 毛 菊 NPP (d) 紫花 针 茅 NPP 
9 * 25 v 2.5 
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图 2 NPP NEP 日 变化 趋势 
Fig.2 Diurnal variation trends of NPP and NEP 
(a) 金 露 梅 NPP (b) 小 葛 草 MPP (c) 风 毛 菊 NPP (d) 紫花 针 茅 NPP 
500 E 300 S 300 ee 
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= 200 = 150 = 150 = 80 
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(e) 青藏 草草 MPP (f) 金 露 梅 NEP (g) 小 项 草 MNEP (bh) 风 毛 菊 NEP 
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(i) 紫花 针 茅 NEP 
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3 NPPNEP 年 变化 趋势 
Fig.3 Annual variation trends of NPP and NEP 
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落 具 有 显著 的 碳 汇 功能 , 且 由 东南 向 西北 碳 汇 功能 
逐渐 减弱 ,其 碳 汇 功能 大 小 排序 为 : 金 露 梅 > 小 葛 草 
> 风 毛 菊 > 紫花 针 茅 > 青藏 蔓草。5 种 植被 VEP 年 空 
间 变 化 趋势 体现 出 ,气温 升 高 和 降水 量 合理 匹配 能 
够 有 效 促 进 植被 的 生长 ,植被 VEP 会 继续 升 高 , 固 
碳 能 力也 会 逐渐 增强 。 

2.3 NPP、NEP 与 气候 因子 的 相关 性 分 析 

2.3.1 NPP NEP 与 降水 的 相关 性 分 析 本 研究 对 
2001 一 2017 年 5 种 植被 的 站 点 年 降雨 量 做 了 统计 ， 
金 露 梅 的 年 平均 降雨 增长 率 最 大 ,但 平均 每 年 不 超 
过 9 mm ,其 余 4 种 群落 的 多 年 年 降雨 趋势 线 大 致 保 
持 一 致 。2001 一 2017 年 三 江 源 5 种 典型 草地 群落 与 
年 降雨 量 相关 性 分 析 得 出 (图 4a ~ e) ,群落 NPP 与 
降水 量 的 相关 性 很 小 , 尼 均 小 于 0.1; NEP 与 降水 量 
的 尽 也 较 小 ,表明 NPP、NEP 受 降水 的 影响 不 显 
著 。NEP 与 降水 的 相关 程度 依次 为 : 风 毛 菊 群 落 
(R=0.123 9) >/] i HEAR (R°=0.04 ) > 28 FEET ae REA 
(R=0.000 6)> 金 露 梅 群落 (R=0.000 3)> 青 藏 蔓草 群 
落 (RR=0.000 03), 值 由 东南 向 西北 逐渐 递减 。5 个 
群落 的 WPP 值 随 着 降水 量 的 增多 有 一 定 的 缓慢 上 
升 趋势 , 且 NEP 也 随 着 降雨 量 的 增多 有 着 细微 的 上 
升 趋势 。 

2.3.2 NPP、NEP 与 气温 的 相关 性 分 析 同 理 对 
2001 一 2017 年 5 种 植被 的 多 年 气温 进行 统计 (图 
5a ~ e) ,三 江 源 5 种 典型 草地 群落 NPP NEP HHA 
温 相关 , 除 金 露 梅 和 风 毛 菊 外 尼 值 略 偏 低 以 外 ,其 


余 3 种 植被 尼 值 均 大 于 0.5 , 故 草 地 群落 与 气温 具有 
显著 相关 性 。 年 平均 气温 为 波动 上 升 的 趋势 ,其 中 青 
藏 蔓草 温度 上 升 最 快 ,每 年 上 升 大 致 0.08%。NPP 与 
气温 的 尼 值 排序 为 紫花 针 茅 群 落 ( 尼 =0.696 8) > 青藏 
蔓草 群落 ( 尼 =0.623 5)> 小 草草 群落 ( 尼 =0.325)> 风 
毛 菊 群落 (R=0.377 5)> 金 露 梅 群落 (RR=0.2163), 具 
有 很 强 的 空间 分 布 特征 。 而 NEP 值 则 与 气温 的 相 
RERA. MPP 与 气温 的 相关 性 整体 上 与 气温 的 
空间 分 布 星相 反 的 趋势 ,其 原因 是 严寒 限制 了 植被 
的 生长 ™, 当 其 限制 条 件 改 善 时 ,其 对 植被 生长 将 
会 具有 较 强 的 促进 作用 。 

NPP. NEP 与 水 热 的 相关 性 对 比 得 出 , NPP、 
NEP 均 受气 温 的 影响 较 大 , 且 成 正 相 关 。 而 降水 与 
它们 的 相关 性 则 不 显著 ,其 规律 与 赵 俊 芳 2 LUYS- 
SAERT™”! 、 裴 志 永 与 等 学 者 的 研究 结果 符合 。 

2.4 TRA FAL (CUE) DH 

碳 利 用 效率 (Carbon use efficiency , CUE) FW] 
级 生产 力 与 总 初级 生产 力 (Gross Primary Productivity, 
GPP) 之 比 ,能 够 反映 植被 吸收 并 固定 CO; 的 能 力 *。 

2017 年 三 江 源 5 个 草地 和 群落 日 碳 利 用 率 均 表现 
为 明显 的 “ 凸 ? 字 形 ( 图 6a ~ e) ,无 明显 波峰 波 谷 。 
空间 上 ,三 江 源 由 东南 向 西北 方向 可 以 看 出 :CUE 
数值 范围 逐渐 由 (-1.0~ 0.8) 缩 小 到 (-0.4 ~ 0.8) , 峰 
值 持续 时 间 大 致 由 240 d 递 减 到 120 d。 表 明 每 年 单 
位 面积 植被 的 固 碳 潜力 也 在 逐渐 下 降 。 

5 种 草地 群落 历年 的 CUE 走势 图 得 出 (图 6f ~ 
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© 400 Ò ooo oo © oroe 
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图 4 NPP、NEP 与 降水 的 相关 性 
Fig. 4 Correlation between NPP(or NEP)and precipitation 
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图 5 NPP WNEP 与 气温 的 相关 性 
Fig. 5 Correlation between NPP (or NEP) and temperature 
(a) 金 露 梅 日 变化 (b) 小 蒿 草 日 变化 (c) 风 毛 菊 日 变化 (d 紫花 针 茅 日 变化 
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图 6 碳 利用 率 日 变化 趋势 及 年 变化 趋势 
Fig.6 Diurnal and annual change trends of CUE 
j, Fa vee GEE CUE WR LIP BA RE 与 降水 呈正 相关 ,与 气温 呈 负 相关 ; 且 随 时 间 的 推移 


针 茅 数值 波动 较 大 ,但 总 体 趋势 是 保持 不 变 ; 小 项 
T MER 、 金 露 梅 皇 明显 下 降 趋 势 。 空 间 上 ,由 东 
南 向 西北 呈现 :缓慢 下 降 - 保 持 稳 定 -缓慢 上 升 的 空 
间 特 征 。 据 大 量 研究 表明 ””: 碳 利用 效率 整体 上 


呈现 逐渐 降低 趋势 ;退化 越 明 显 的 地 区 CUE 的 值 越 
小 。 东 南部 3 个 草地 群落 CUE 变化 趋势 同 现 有 ” 研 
究 一 致 ,主要 与 气温 呈 负 相关 , 主要 原因 是 5 种 植被 
多 年 降雨 趋势 基本 不 变 , 气 温 都 表现 为 明显 增加 趋 
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势 , 那 么 高 寒 地 区 的 植被 在 气温 升 高 并 不 显著 的 情 
况 下 ,植被 的 生长 仍然 受到 限制 ,因此 CUE 的 变化 
趋势 不 明显 。 


3 结论 


(1) 三 江 源 典型 草地 群落 具有 碳 汇 作用 , 且 
2001 一 2017 年 间 碳 汇 功能 在 增强 。 三 江 源 的 5 种 草 
地 植被 一 年 中 的 活跃 生长 期 在 第 50 ~ 300 d 这 个 时 
段 ;各 个 群落 的 NPP 值 近 17 a 来 均 表 现 为 上 升 趋 
势 , 其 中 青藏 蔓草 WPP 值 最 低 , 为 55.93gC.m?.a ， 
相 较 之 下 不 及 金 露 梅 NPP 值 (342.00 g C.m-.a-0) 
的 十 分 之 一 ;整个 研究 区 年 平均 NEP 值 均 在 20 gC 
ma , 均 大 于 0。 

(2) BEY CUE MEIJER ,小 草草 .紫花 针 茅 、 风 
毛 菊 、 青 藏 蔓草 CUE 值 均 在 0.625 以 上 ,有 较 强 固 碳 
潜力 ; 且 空 间 上 东南 部 区 域 植被 固 碳 潜力 逐渐 降 
低 , 西 北部 逐渐 增加 。 

(3) NPP 与 气温 呈 显 著 正 相关 。 其 中 紫花 针 茅 
群落 (R=0.696 8) 青藏 草草 群落 (及 =0.623 5) 和 小 
ERPS (R°=0.525) PB RIYKFO.S. 

(4) 本 人 研究 中 很 多 数据 均 体 现 出 独特 的 空间 分 
布 格局 。 由 东南 向 西北 纵 观 三 江 源 :高 程 由 低 到 
高 ,水 热 组 合 逐 渐 减 少 ,因此 植被 由 高 大 密集 的 灌 
木 . 乔 木 向 稀 琉 矮小 的 耐寒 耐 旱 高 山 植 被 过 渡 , 植 
W NPP VEP 数 值 逐 渐 递 减 , 固 碳 潜力 逐渐 增强 。 
另外 ,本 论文 中 BIOME-BGC 模型 的 参数 优化 设置 
有 待 深入 研究 。 
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Simulation and dynamic change of carbon source/sink in the typical grassland 
communities in the Three River Source Area from 2001 to 2017 


CHEN Xue-jiao, ZHOU Wei, YANG Han 
(Chongqing Jiaotong University, College of Architecture and Urban Planning, Chongqing 400074, China) 


Abstract: In this paper, five typical grassland communities (Potentilla fruticosa, Stipa purpurea, Saussurea 
japonica, Kobresia pygmaea, and Carex moorcroftii communities) in the Three River Source Area, Qinghai Province, 
China were selected as research objects. Basic geographic data, meteorological data, and vegetation physiological 
parameters data were used to estimate the net primary productivity (VPP) and the net ecosystem productivity (VEP) 
in these grassland communities from 2001 to 2017 based on the BIOME-BGC model. This paper also investigated 
the change characteristics of NPP and NEP and their correlation with temperature and precipitation. The 
characteristics of the change of the carbon use efficiency (CUE) were also explored. Results showed that the spatial 
pattern of the values of NPP and NEP in the Three River Source Area reduced gradually from southeast to 
northwest. The multi-year average value of NPP in 5 typical grassland communities was 196.06 g C +m?» a~, which 
showed an increasing trend year by year. Among them, the annual average value of NPP in Potentilla fruticosa was 
the highest (342.00 g C +m” a '), and the annual average value of NPP in Carex moorcroftii seagrass was the lowest 
(55.93 g C+-m a). The annual mean value of NEP of the five grassland communities was 49.02 g C+m™’+a'. The 
NEP values of Potentilla fruticosa, Stipa purpurea and Carex moorcroftii increased slowly, while the NEP values of 
Saussurea japonica and Kobresia pygmaea decreased slowly. Therefore, the grassland ecosystem in the Three River 
Source Area has significant carbon sequestration function of absorbing CO2. A significant positive correlation 
between NPP and temperature was observed, while a low correlation between NPP, NEP, and precipitation existed. 
All the five communities have great potential for carbon sequestration. The carbon utilization rate of other vegetation 
communities was above 0.625, except the Potentilla fruticosa. 


Key words: carbon source/sink; NPP; NEP; BIOME-BGC model; climate factor; Three River Source Area 


